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Observational consequences of the tidal disruption of stars (TDEs) by supermassive black holes 
(SMBHs) can enable us to discover quiescent SMBHs and constrain their mass function. More¬ 
over, observing jetted TDEs (from previously non-active galaxies) provides us with a new means 
of studying the early phases of jet formation and evolution in an otherwise “pristine” environ¬ 
ment. Although several (tens) TDEs have been discovered since 1999, only two jetted TDEs have 
been recently discovered in hard X-rays, and only one. Swift J1644+57, has a precise localization 
which further supports the TDE interpretation. These events alone are not sufficient to address 
those science issues, which require a substantial increase of the current sample. Despite the way 
they were discovered, the highest discovery potential for jetted TDEs is not held by current and 
up-coming X-ray instmments, which will yield only a few to a few tens events per year. In fact, 
the best strategy is to use the Square Kilometer Array to detect TDEs and trigger multi-wavelength 
follow-ups, yielding hundreds candidates per year, up to z ~ 2. Radio and X-ray synergy, how¬ 
ever, can in principle constrain important quantities such as the absolute rate of jetted TDEs, their 
jet power, bulk Lorentz factor, the black hole mass function, and perhaps discover massive black 
holes (MBH) with < 10^Mq. Finally, when comparing SKA results with information from optical 
surveys like LSST, one can more directly constrain the efficiency of jet production. 
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1. Introduction 

Mounting observational evidence is supporting a scenario where most galactic nuclei host 
supermassive black holes (SMBHs). Gas inflow from larger scales causes a small (~ 1%) fraction 
of SMBHs to accrete continuously for millions of years and shine as quasars. Most of them are 
instead “quiescent”, accreting -if at all- at a highly sub-Eddington rate. Observationally, it is 
therefore hard to assess the presence and measure the mass of most of SMBHs, beyond the local 
universe. Occasionally, however, a sudden increase of the accretion rate may occur, if a large mass 
of gas, e.g. a star, falls into the tidal sphere of influence of the black hole and And its-self torn apart 
and accreted. We call these events “tidal disruption events” (TDEs). TDEs can result in sudden 
flares: these can reach the luminosity of a quasar but they are rare (~ 10“^yr“^ per galaxy) and last 
for only a few up to several months. We will show in this Chapter how the SKA sensitivity and its 
large held of view could increase the TDE sample and have a large discovery potential. 

The detection and study of these flares can deliver other important astrophysical information, 
beyond flagging the presence of a SMBH. The theoretical expectation is that TDE flares would 
be caused by sudden accretion of the star debris. If the star is completely disrupted, its debris is 
accreted at a decreasing rate of M (Rees 1988; Phinney 1989). Therefore, TDEs allow 

us to study the formation of a transient accretion disc and its continuous transition through dif¬ 
ferent accretion states. As the accretion rate decreases, we can in principle observe a disc in an 
initial super-Eddington phase, lasting several months; the disc then would become slim and later 
enter in a thin disc regime. Anally ending its life, years later, in a radiative inefficient state. The 
super-Eddington phase — which occurs only for SMBH masses M 10^ M© — is theoretically 
uncertain, but it may be associated with a powerful radiative driven wind (Rossi & Begelman 2009), 
which thermally emits ~ 10^^ — 10^^ erg s“^, mainly at optical frequencies (Strubbe & Quataert 
2009; Lodato & Rossi 2011). The disc luminosity (~ 10^^ — 10“^^ erg s“^), instead, peaks in the 
far-UV/soft X-rays (Lodato & Rossi 2011). Of great theoretical importance would also be the 
possibility to observe the formation and evolution of an associated jet, powered by this sudden 
accretion episode. 

The first TDEs were discovered in ROSAT surveys of the X-ray sky (e.g., Komossa & Bade 
1999; Grupe et al. 1999) — see Komossa (2002) for a review. Later, GALEX allowed for the 
selection of TDE at UV frequencies (Gezari et al. 2009, 2012); many of the most recent TDE 
candidates are now found in optical transient surveys (van Velzen et al. 2011a; Cenko et al. 2012; 
Chomock et al. 2014; Arcavi et al. 2014). An alternative method to select TDE candidate is to 
look for optical spectra with extreme coronal lines (Komossa et al. 2008; Wang et al. 2012). These 
“thermal" features are believed to be associated with emission from the disc or the radiative driven 
wind. 

Moreover, recently the Swift Burst Alert Telescope (BAT) instrument (Burrows et al. 2011; 
Cenko et al. 2012), triggered two TDE candidates in the hard X-ray band. A multi-frequency 
follow-up from radio to y-rays revealed a new class of non-thermal TDEs. It is widely believed 
that emission from a relativistic jet is responsible for the hard X-ray spectrum (with power-law 
photon index /3 ~ 1.6 — 1.8) and the increasing radio activity (Levan et al. 2011), detected a few 
days after the trigger. The best studied of the two events is Swift J1644-1-57 (Sw J1644 in short), 
whose radio lightcurve (see next section) is used in this Chapter to estimate detection rates in the 
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radio band. The two main features that support the claim that Sw J1644 is a tidal disruption event 
are i) the X-ray lightcurve behaviour, that follows after a few days from the trigger, and 

ii) the radio localization of the event within 150 pc from the centre of a known quiescent galactic 
nucleus (Zauderer et al. 2011). 

1.1 Radio emission of jetted TDEs 

In 2011, BAT discovered a peculiar gamma-ray burst, now known as Sw J1644. What first 
caught the attention was the flaring activity, which caused multiple BAT triggers, followed by an 
X-ray light curve that extended to much longer timescales than any known GRB (Levan et al. 2011). 
The source was conhrmed to be of extragalactic origin, its peak luminosity ~ 100 times higher than 
bright active galactic nuclei (Burrows et al. 2011; Levan et al. 2011). Follow-up radio observations 
revealed a transient radio counterpart (Zauderer et al. 2011) and the multi-wavelength SED of the 
source resembled that of BL Lac objects (Bloom et al. 2011). Key evidence that Sw J1644 is, in 
fact, a newly launched jet, is the rapid X-ray variability implying Doppler-beamed emitting regions 
(Bloom et al. 2011). A second BAT-triggered TDE candidate. Swift 2058-1-05, that also has a 
transient radio counterpart was discovered by Cenko et al. (2012). 

A rich dataset of radio follow-up observations of Sw J1644 is presented in Berger et al. (2012); 
Zauderer et al. (2013). At 5 GHz, the peak in the light curve occurs about 100 days after the hrst 
gamma-ray trigger, at Fy = 15 mJy (vLy = 4 x 10“^^ergs) (cf. Eig. 3 in Zauderer et al. 2013). 
During the hrst years the self-absorption frequency is around 10 GHz; at 2 GHz the peak of the 
light curve is expected to occurs about 500 days after the stellar disruption. The lightcurve at 1.4 
GHz is shown in Eig. 3.2. Data were collected from 5 days after the Swift trigger, with a hux 
of Fy > 0.2 mJy. The radio light curve of the second relativistic TDE, Swift 2058-1-05, has been 
sampled for only a few epochs (at 5 GHz the hux is 0.9 mJy). Its luminosity is similar to Sw J1644, 
but the spectrum is rather hat (Fy V®) from 4.5 GHz up to 22.5 GHz Cenko et al. (2012). 

The light curve of Sw J1644 can be modeled as synchrotron emission from the forward shock 
(Giannios & Metzger 2011; Metzger, Giannios & Mimica 2012), although this model requires an 
increase of the isotropic-equivalent jet kinetic energy (Berger et al. 2012), which is much stronger 
than what is expected from the fallback rate of stellar debris (see also Kumar et al. 2013). Interest¬ 
ingly, the X-ray light curve Sw J1644, shows a very steep drop when the (theoretically estimated) 
fallback rate crosses the Eddington limit (Zauderer et al. 2013). This could imply that the jet engine 
has switched off due to a state change of the accretion disk, as suggested by the hrst models of TDE 
jets (Giannios & Metzger 2011; van Velzen, Kdrding & Ealcke 2011). 

2. Constraints on TDE jets from past and present radio surveys 

So far, surveys of radio transients have yielded a few events (e.g. Gal-Yam et al. 2006; Ban¬ 
nister et al. 2011) — see Erail et al. (2012) for a review. However, due to the lack of follow-up 
observations it has been difficult to test if these are TDE jets. Eor a given light curve model of the 
radio emission from TDE jets, the upper limits on the event rate from these surveys can be used to 
constrain the rate of these events. 

As shown in van Velzen et al. (2013), the observed light curve of Sw J1644 can be used to 
make an order of magnitude estimate of snapshot rate of events that have a similar light curve. 
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Input to this estimate are the per-galaxy rate of tidal disruption jets (A^tde) and the mean volume 
density of black holes Pbh that create observable TDE jets. If we assume that only jets with an 
inclination smaller than 1 /Fj are visible, we arrive at the following expression for the snapshot rate 
above a given flux limit 


R(Fv. 


. 1 r\—3 1 


\ 3/2 A rj-,-\T 


(2.1) 


Here AT denotes the time interval along which Sw J1644 flux is above Here we normalized 
snapshot rate assuming that the rate of jetted TDE iVrDJ is similar to the rate of thermal TDE 
(we return to this assumption below). The rate of thermal TDE can be estimated from soft X- 
ray (Donley et al. 2002) or optical (van Velzen & Earrar 2014) surveys, a conservative value is 
^TDE ~ 10 ^ yr '. 

Eor AT = 1 yr, the light curve of Sw J1644 implies Fy^sw = 10 mJy, 2 mJy at 5 and 1 GHz, re¬ 
spectively. Adopting r,- = 2 as estimated from modelling the late-time radio emission of Sw J1644 
(Berger et al. 2012), we obtain our estimate of snapshot rate as shown in Eig. 1. This rate is close 
to the existing upper limits of radio transient surveys. 

In Eig. 1 we also show the predicted rate at 1 GHz based on the model of van Velzen, Kording 
& Ealcke (2011). In the optimistic scenario of these authors, 20% of the accretion power is con¬ 
verted into jet power (i.e., in analogy with AGNs, tidal disruption jets are always radio-loud). The 
synchrotron luminosity is computed using the conical equipartition jet model (Blandford & Konigl 
1979; Ealcke & Biermann 1995), which yields a peak luminosity (for on-axis events) that is lower 
than the peak of Sw J1644 (hence the predicted snapshot rate is also lower). In the conservative 
scenario of van Velzen, Kording & Ealcke (2011), tidal disruption jets are only radio-loud when 
the fallback rate drops below 2% of the Eddington limit. 


3, SKA as the best hunter of jetted TDEs 

In this section we report on the rate of TDEs expected to be detected with SKA operating 
in survey mode. The SKA 1-survey (SKAl-SUR) provides us with an optimal tool to catch and 
discover jetted TDEs. We adopted a strategy which relies on a trade-off between sky coverage 
and sensitivity, fully compliant with the SKAl-SUR science requirements reported in the Baseline 
Design Document (Dewdney et al. 2013). We request to achieve an half sky coverage (20,000 
deg^) with a 2-day cadence at a 5 — (7 flux limit of 90 fjJy. This can be achieved by SKAl-SUR. 
SKAl-SUR can cover this area with 1111 pointings at 1.4 GHz (the SKAl-SUR field of view at 
this frequency is ISdeg^) at a sensitivity of 3.72 microJy hr“^/^ (see Table 1 in Dewdney et al. 
2013), or at a flux limit of 90 fjJy in 2 days. 

3.1 Snapshot rate estimate 

We first make a rough estimate based on Eq. 2.1. At 1 GHz, the observed light curve of 
Sw J1644 is above 2 mJy for less than one year (AT ss O.Syr) (Zauderer et al. 2013), for a flux limit 
90 jjJy we thus obtain a snapshot rate of 0.2 deg“^ (Fig. 1), implying w 4000 events for the SKA 
transient survey adopted here if Fj = 2 is assumed. 
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■ ■ van Velzen et al. (2011): Optimistic (1 GHz) 
. van Velzen et al. (2011): Conservative (1 GHz) 
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Figure 1: Predicted snapshot rate and upper limits from current and future radio surveys. The predicted rate 
for Sw J1644 is based on eq. 2.1 and assumes a 100% jet production efficiency, i.e., every TDE makes a 
jet with a light curve similar to Sw J1644. The blue, red, and magenta triangles show (potential) 2-o upper 
limits at 5, 3 and 1 GHz, respectively. References: B7/F12: Bower et al. (2007); Frail et al. (2012); de Vries 
et al. (2004), FIRST-NVSS: Gal-Yam et al. (2006); Croft et al. (2010); Scott (1996); Mooley et al. (2014). 

3.2 Radio Lightcurve and Monte Carlo simulations 

We here extend our rate calculations to explicitly account for redshift, BH mass M and stel¬ 
lar mass m* dependences, using Monte Carlo simulations (see Donnarumma & Rossi 2014, for 
details). The radio lightcurve at 1.4 GHz of Sw J1644 is plotted in Fig. 3.2. To investigate the 
detection capability of SKAl-SUR, we assume that all events have the same light curve shape as 
Sw J1644. We model the mission as synchrotron radiation and use the peak flux and characteristic 
frequencies measured at different times by Berger et al. (2012, see their eq.4 and dash-dotted line 
in our Fig.3.2) extended at late epochs by Zauderer et al. (2013). The lightcurve behaviour in the 
first ~ 8 days (where there are no data) has been extrapolated from later epochs. The result is a 
quite sharp rise, that can be considered a conservative assumption: a more gentle rise would ensure 
more detections. Note that ours is not a fit to the data, but it is a modelling. From Fig. 2 it is clear 
that the radio flux does not show the temporal characteristic behaviour oc of a tidal disruption 
event, that it is present in the X-ray band. Therefore, a follow up strategy at higher frequency is 
required for identification (see Section 3.5). 

An appropriate spectral correction can describe events which explode at different redshifts 
than Sw J1644 ^ Rescaling the flux for different BH and stellar masses requires instead to assume 
a jet model. A first possibility is to describe the jet evolution with a Blandford McKee solution 

* We here assume that all jets have approximately the same Doppler factor 5 ~ T^- ~ 2, as measured in radio for Sw 
J1644 


5 





Jetted TDEs 


I. Donnarumma 





10.000 


1.000 


0.100 


0.010 


0.001 



100 200 300 400 

t - T, (days) 


Figure 2: The radio (1.4 GHz) light curve of Sw J1644: data (black circles) from Berger et al. (2012); Za- 
uderer et al. (2013), modelling (dash-dotted black line) from Donnarumma & Rossi (2014). Figure adapted 
from Donnarumma & Rossi (2014). Note that the black line is not a fit. We trace it by assuming synchrotron 
emission and the characteristic frequencies and peak flux measured by Berger et al. (2012); Zauderer et 
al. (2013) at different times. The result slightly underpredicts the data, but recovers the overall temporal 
behaviour which is rather flat. 


(Blandford & Mckee 1977), usually adopted for y-ray burst afterglows. In this case, the radio flux 
does not explicitly depend on the BH mass. A second approach is to assume that the luminosity 
is proportional to the peak of the jet power, which in turn is proportional to the peak fall-back rate 
Lj p oc Mp. This introduces a mass dependence in the luminosity as Ly Lj.p . Here¬ 

after we will focus on this second BH mass dependent lumimosity (MDL) model and we refer the 
readers to Donnarumma & Rossi (2014) for a comparison with the Blandford Mckee solution. The 
stellar mass distribution in galactic nuclei may vary, but there is no solid observational evidence. 
We therefore assume a ”universal“ Kroupa initial mass function Kroupa (2001), as observed in our 
Galaxy. The BH mass function is instead explained more in detail in the following section. 

3.3 Black hole mass functions 

The mass distribution of black holes as a function of redshift is a key ingredient to calculate 
TDE rates. We adopt the BH mass function (MF) predictions according to the models labeled G and 
G(z) in Shankar, Weinberg & Miralda-Escude (2013). In practice, we consider the two accretion 
models which yield the largest and the lowest MF predictions in order to reflect the uncertainty 
due to the black hole mass distribution on our expected TDE rates. The “intrinsic” TDE rate as a 
function of redshift is defined by 

rttimax 

R{z)= (^{M,z)NrDEV{z)dM, (3.1) 

where Atde = 10“^ yr“^ is as our fiducial intrinsic TDE rale, ^{M,z) is Ihe black hole comoving 
number density and V (z) Ihe comoving cosmological volume. In our calculalions, we consider 
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only BH masses above M^in = 10® M©, because MFs are not observationally constrained below 
that mass. 

3.4 Monte Carlo simulations 

Assuming the radio modellings described in section 3.2, we performed Monte Carlo simula¬ 
tions (MCs) to derive the rate of jetted TDEs expected to be detected by SKAl-SUR operating in 
wide survey mode . 

The intrinsic rate R{z) (eq. 3.1) is reduced by a factor of 1 /2 T? for the relativistic beaming. 
Our fiducial value is T, = 2 (Zauderer et al. 2013; Berger et al. 2012), but our results may be 
easily rescaled by assuming different values of T^. Finally, a further scaling factor accounts for the 
fraction of the sky accessible at each time by the instrument. As already discussed in section 3.3, 
we will provide a range of rates whose limits reflect the uncertainty in the BH ME predictions. 

The specific observafional sfrafegy for SKAl-SUR explained al Ihe beginning of Ihis seclion 
assumes an half sky coverage (20,000 deg^) wilh a 2-day cadence al a 5 — C7 flux limil of 90 fjJy. 

For each evenl, wilh given redshifl, BH and stellar mass, Ihe lighl curve is predicted as de¬ 
scribed in seclion 3.1. We exlracl Ihe Irigger lime randomly from a 1 year uniform dislribulion, 
and we calculate Ihe average flux over 2 days since Ihe Irigger and Ihen we compare il wilh Ihe 
SKAl-SUR flux limil fixed for our slralegy. 

We show in Fig. 3, Ihe dislribulion of Ihe TDE rate as a funclion of z (righl panel) and Ihe 
corresponding BH mass dislribulion (left panel) for Ihe Iwo BH MFs assumed (solid and doited 
blue lines). In Ihe righl panel, we overplol Ihe rate dislribulions for TDEs wilh BH masses lower 
lhan 10^ Mq (solid and doited red lines). The lolal rate dislribulion peaks around z ~ 0.5 and il is 
dominated by lighl SMBHs locally and for z above Ihe peak. The maximum redshifl lhal SKA can 
allain is close lo z ~ 2. 

The lolal number of objecls oblained by inlegraling Ihese dislribulions in z and M^^^ is ~ 300 
and ~ 800 yr“* for Ihe G and G(z) models, respeclively. We note lhal our integrated rates are 
signihcanlly lower lhan Ihe hrsl simple eslimales reported in seclion 3. This is mainly due lo 
Ihe fad lhal we accounted for Ihe stellar mass dependence, lhal inlroduces a large varialion in Ihe 
inlrinsic luminosity of Ihese evenls. Neverlheless, our resulls imply lhal a signihcanl increase in Ihe 
number of detected jelled TDEs can be achieved lhanks lo our slralegy in Ihe framework of SKAl- 
SUR. However, as menlioned before, Ihe TDE signalure seems lo appear al higher frequencies. We 
Iherefore proceed lo consider X-ray facililies. 

3.5 Radio and X-ray synergy 

Fulure X-ray surveys can deled much fewer candidates lhan SKA. The e-Rosila ah sky survey 
(Merloni el al. 2012) is expected lo Irigger a few jelled TDEs by using ils “hard" X-ray band (2 — 10 
keV). Calculalions based on Ihe Swift slalislics (by using SW 2058 as a prototype) predicted 1 
objecl to be detected per scan (6-monlh long) probing Ihe universe up to z ~ 4.5 (Khabibuhin el al. 
2014). 

We also envisage Ihe use of X-ray telescopes to follow up radio candidates and conhrm Iheir 
nalure. To quanlilalively explore Ihis possibility, we model Ihe X-ray emission of Sw J1644 in 
Ihe 1-10 keV band, assuming lhal Ihe X-ray luminosity scales as Ihe fall back rate. Lx M. (see 
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Figure 3: Rate of events detected by SKA as a function of redshift (right panel) and their mass distribution 
(expressed in percentage, left panel), for two different black hole mass distribution functions (solid and 
dotted blue lines). The rates for BH masses lower than 10^ M© are also shown (red lines, right panel). 
Figure adapted from Donnarumma & Rossi (2014). 


Donnarumma & Rossi 2014, for a full description of the methodology). In our estimates, we 
assumed that a follow-up can be arranged within 1 day from the radio trigger. To claim detection, 
we require that the observed flux can be detected with a signal to noise of at least 10 during four 
pointings, spread over four days. This requirement implies that a flux threshold (unabsorbed flux) 
of ~ 5 X 10“^^ erg cm“^ s“* in the 1-10 keV band is in principle enough to identify all the SKA 
selected candidates. 

Therefore, given their sensitivity, future X-ray experiments like Athena (Nandra et al. 2013) 
and a LOFT-like mission (Feroci et al. 2012) could offer a unique chance to follow-up and charac¬ 
terize SKA-triggered TDEs. If Athena will have a follow-up efficiency of 50% associated with the 
sky accessibility, it will be able to identify about 50% of all the radio sample. For a LOFT-like mis¬ 
sion, we would expect that a sub-sample of radio TDEs between ~ 100 and ~ 350 will be detected 
per year up to z ~ 1.2 by observing with its main instrument, the Large Area Detector. 


4, Constraints on the jet production efficiency 

As recalled above, only two jetted TDEs have been detected, while detections of the thermal 
candidates, related to the presence of an accretion disk, have been more numerous. The question 
then arises whether this is due to observational biases or to a small jet opening angle or to an 
intrinsic low efficiency of transient accretion disks to produce (luminous) jets. 
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To answer this question, one may first compare predictions in the Swift/BAT bands with the 
observed rate of ~ 0.3 yr“*. It turns out that the most pessimistic predictions by Donnarumma & 
Rossi (2014) are still a factor of tens higher than observed and it would be tempting to reconcile 
this discrepancy by invoking a jet production efficiency of a few percent. However, it has been 
shown that this inference cannot be confidently drawn (Donnarumma & Rossi 2014). On the one 
hand, other uncertain parameters -like the value of the beaming factor in the X-ray emitting region 
- can mitigate the discrepancy. On the other, the modelling of the characteristics (e.g. flux limit 
and sky coverage) of an effective Swift/BAT survey (in fact. Swift is not operating in survey mode) 
it is not at all straightforward and does not allow for firm predicfions. Plus, with just two events, 
the Poisson uncertainty is very large. 

Another approach for estimating the jet production efficiency is fo use instead the observed 
thermal candidates. Thermal TDEs may also contain jets. In fact, based on an analogy with X- 
ray binaries, one may anticipate that each TDE is accompanied by a jet (see van Velzen, Kording 
& Ealcke 2011, and references therein). Because they are observed off-axis, these jets are less 
luminous and will not trigger Swift, but they can be detected through radio or X-ray follow-up 
observations. 

The existing radio follow-up observations of known thermal TDEs are summarized in van 
Velzen et al. (2013). So far, none of the optical/UV-selected TDE have been detected at radio 
frequencies (see also Arcavi et al. 2014). Using a simple model to scale the observed light curve 
of Sw J1644 for off-axis observed we find that the current radio observations are deep enough to 
detect jets similar to Sw J1644 to an inclination of i < 70° (van Velzen et al. 2013). We thus know 
that the jet production efficiency has to be less than 100%. 

Two of the seven TDE candidates that were detected by ROSAT or XMM have been detected 
at 5 GHz in follow-up observation by Bower et al. (2013), namely IC 3599 (Komossa & Bade 
1999) and RX J1420-I-5334 (Greiner et al. 2000). The radio observations were obtained about ten 
years after the estimated time of disruption, their transient nature remains to be confirmed (in fact 
one of the radio-detected sources, IC 3599, is more likely to be an AGN). 

We can thus conclude that none of the thermal TDE show strong evidence for relativistic jets. 
This could imply that the fraction of TDE that launch jets similar to Sw J1644 less than « 10%, 
although currently is not clear if the detection of optical/UV emission implies a bias against the 
detection of radio emission. To be able to detect optical or UV emission, one requires a circum- 
nuclear environment with little absorption by dust, while the radio luminosity can be boosted by a 
high density circumnuclear gas. 

Present constraints on the jet production efficiency suffer of the poor statistics due to the 
small X-ray and radio samples. The radio rates of jetted TDEs predicted in this work, will help to 
overcome the present limitation due to the very large sample expected to be triggered by the SKA. 
Besides serving as a trigger, the high survey speed of the SKA implies that it will be possible to 
have near-simultaneous observations with LSST surveys. The detection rate of optical TDE with 
LSST is ~ 1000yr“^ (e.g., Gezari et al. 2009; van Velzen et al. 2011a), hence the SKA observations 
will be able to measure a jet production efficiency for thermal TDE as small as 0.1%. Comparing 
the SKA-triggered TDEs with the SKA observations of thermal TDE discovered by LSST will 
provide a new window to study the jet physics of massive black holes. 
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5. Implications for the cosmological growth of SMBHs 

The observed properties of central SMBH and their hosts show evidence for a strong link be¬ 
tween SMBH and galaxy formation (Ferrarese & Merritt 2000; Gebhardt et al. 2000). Outstanding 
questions include how the massive black hole seeds (MBH, ~ 10^ — 10^ M©) were initially formed, 
and their successive mass growth, which determine their mass functions (Volonteri 2010). TDEs 
will provide a great opportunity in this quest. 

Fig.3 shows that SKA will more easily detect light (< lO^M©) MBHs locally and at z > 1. 
Since the lighter the more luminous, this result suggests the possibility to observe BHs in the 
elusive "intermediate mass regime" in the distant Universe. 

As described in section 3, in this Chapter we focused our work on BHs with masses > lO^M©. 
This because of the very high uncertainties affecting the lower mass side of the BH mass functions. 
However, the recently reported TDE candidate from a 10^-^ Mq black hole (Donato et al. 2014) 
shows the great potential to exploit TDEs to seek out MBHs in the mass range of 10^ - 10^ Mq, 
with important consequences for BH formation models. These black holes must exist in low- 
mass dwarf galaxies or perhaps as remnants inside massive galaxies. Alternatively, they may have 
formed in the centre of globular clusters^ (Maccarone 2004). 

We would like to point out, that the high-precision location of the radio-TDE with respect to 
the galaxy core will be of great relevance in the understanding of both SMBH and IMBH formation. 
TDEs might reveal recoiling black holes (e.g. Komossa 2012) that have wandered slightly away 
from the nucleus, off-nuclear SMBH due to ongoing merger activity (Comerford & Greene 2014), 
or for example black holes that were formed in the centre of globular clusters (Maccarone 2004). 
In the case of galaxy mergers, the stellar tidal disruption rate can be temporarily strongly enhanced 
(e.g., Chen et al. 2009). Such TDE would point us to an ongoing merger, but would likely be off-set 
from the geometrical center of the merging system. In the case of recoiling BHs due to e.g. SMBH 
coalescence, gravitational wave recoil will insert a kick velocity on the newly formed single black 
hole, which will carry a retinue of stars with it, as it is ejected from the core of its host galaxy. 
These bound stars will undergo tidal disruption (Komossa & Merritt 2008), leading to a radio event 
offset from the galaxy core. 

6. Observational strategies 

As showed in Section 3, SKA can play an important role in discovering jetted TDEs. To 
achieve this, it will be necessary to find TDEs by means of a commensal search in all SKA ob¬ 
servations, or by carrying out a near real-time transient search in a high cadence dedicated survey 
observations. The latter is preferred since our goal is to detect TDEs early on in order to allow for 
the necessary follow-ups (to confirm fhe TDE origin). The real advanfage of fhis sfrafegy is fhaf if 
will be able fo build up a large sample of jelled TDEs up lo z ~ 2, where fhe low-mass populalion 
of SMBHs could dominate, providing a unique opporlunily lo conslrain fhe low mass end of fhe 
SMBH mass funclion in fhe far Universe. 

^Off-nuclear BH in the mass range ~ 10^ — 10^ Mq, as a likely power source of the most extreme ultra-luminous 
X-ray sources, and the ones that are believed to have formed in the centre of globular clusters are often referred to as 
intermediate-mass black holes (IMBH) in the literature. 
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The SKAl-SUR will provide us with an optimal instrument to detect transient events as TDEs. 
Our strategy, explained in Section 3, is compliant with the Baseline Design Document (Dewdney et 
al. 2013) and predicts a few hundred detections with X-ray follow up for identification. However, 
a significant sample could also be collected, assuming a 50% reduction in sensitivity due to re¬ 
baselining. This is because the 1.4 GHz light curve of Sw J1644 increases with time at least up to 
600 days (Zauderer et al. 2013), therefore a multiple cadence survey with longer integration times 
(~ 8 days) needed to achieve the same flux limit of our strategy should not result in a significant 
loss of radio detections. However, radio triggers at late times could be problematic for follow-up 
observations, because sources would get fainter in both optical and X-ray energy range, preventing 
the identification of higher redshift TDEs. Nevertheless, the bulk of the sample at z < 0.5 should 
still be successfully followed up at higher frequencies. 

On the basis of these considerations, TDE can be commensally searched for, by exploiting 
the SKAl-SUR all-sky survey (31,000 deg^ at 2 ^uJy beam“' at 1.5-2 arcsec resolution) proposed 
in the framework of both the Magnetism and Continuum science cases (see Johnston-Hollit et al. 
2015; Prandoni et al. 2015). If carried out by visiting individual fields mulfiple limes wilh a 
cadence of ~ 10 days over Ihe Iwo year period, if will be very efficienl in finding exlra-galaclic 
Iransienls down lo ~90 ^uJy beam“ ^ . 

As menlioned before, Ihe 1.4 GHz lighl curve of our prololype TDE does nol bear any specific 
foolprinls of a TDE nalure. Thus we cannol use only ils properties or variability lo dislinguish a 
TDE candidate from neilher a slowly variable AGN or a GRB. However, we can accounl for AGN 
conlaminalion by cross-correlaling Ihe Iransienl posilions wilh deep AGN calalogues. A significanl 
advance in Ibis direclion will be achieved by Ihe large amounl of dala lhal LSST will provide in 
Ihe near fulure logelher wilh Ihose collecled after one year of SKA and SKA precursor (ASKAP) 
observalions. On Ihe olher hand, precise localizalion of Ihe radio Iransienl in Ihe core of galaclic 
nuclei will help assessing Ihe nuclear origin: Ihe belter Ihe localisalion, Ihe slronger is Ihe case for 
a TDE. 

This means lhal a rapid oplical follow-up aimed al determining Ihe hosl galaxy has lo be 
foreseen. Once Ihe hosl galaxy is found, lhanks lo a resolulion of aboul 2 arcsec (SKAl-SUR) and 
0.6 arcsec or belter (SKAl-MID) (Dewdney el al. 2013) il will be possible lo localize Ihe brighter 
Iransienls wilh a precision of ~ 100 milliarcsecond (mas), essenlial lo separate nuclear Iransienls 
from olher phenomena (e.g., GRB). This angular localizalion is al Ihe currenl alignmenl accuracy 
level of radio and oplical reference frames (Orosz & Erey 2013), corresponding lo linear scales of 
0.5-0.8 kpc al z > 0.3. 

We note lhal by delecling radio emission associated wilh slar formalion, pre-flare SKA obser¬ 
valions can help lo idenlify jelled TDEs. By combining Ihe individual snapshol oblained after ~ 1 
year of observalions, we expecl lo be able lo see typical slar-forming galaxies (L lO^O W/Hz; 
Condon (1992)) up lo z ~ 0.5. Al Ihis redshift we also expecl lo see mosl of our detected jelled 
TDEs (see Eig. 2), hence mosl TDEs lhal occur in slar-forming galaxies will have a detected hosl 
galaxy. These hosl deleclions will allow for an accurate localizalion of Ihe radio Iransienl wilh re- 
specl lo Ihe center of ils hosl galaxy, which will help lo rejecl SNe and GRBs (TDE from SMBHs 
occur in Ihe nucleus of Ihe galaxy, while SNe/GBRs Irack Ihe stellar density). The idenlilicalion of 
nuclear radio Iransienls will allow for faster Iriggers for follow-ups in X-rays (by avoiding Ihe need 
for Ihe oplical follow-up for Ihe hosl idenlilicalion), which will give us a firm signalure of Ihe TDE 
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nature if the trigger occur within 1 year from the beginning of the event (see section 3). A synergy 
with SKA and VLBI will help to characterize the jet formation in the pristine environment of the 
TDE hosts (Paragi et al. 2014b). 

Given the huge amount of TDE candidates predicted and the expensive strategy, we do not 
need to know immediately for all hundreds of events, what they are. In the first stage, we could 
concentrate the efforts on the brightest events (tens of objects) and for them activate follow-up 
observations at higher energies. The rest of those events could be analyzed later by means of 
commensal search in SKA observations. 

Einally, SKA will play an important role even for the follow-up of all TDEs detected by the 
Gaia mission (expected to end soon before SKA 1-MID early operations) thanks to the SKAl- 
MID sub-arcsec angular resolution. Immediate follow-up observations of TDEs detected in the 
X-rays and by LSST in the optical will be important to constrain the jet production efficiency af 
low redshifl. This will require flexible operafions and fhe abilify fo respond fo friggers from ofher 
insfrumenfs (see Tender 2015). 

7. Conclusions 

In fhis Chapfer, we have shown how SKA alone, fhanks fo ifs combinafion of survey speed 
(SKAl-SUR, for fransienf searches) and sub-arcsec angular resolufion (SKA 1-MID, for TDE iden- 
lificalions), and in synergy wifh coeval insfrumenfs af higher energies, can play a leading role in 
discovering jelled TDEs and allowing for fhe firsl lime slalislical sludies of Iheir fealures. We pre- 
dicl a sample of hundreds, up fo redshifl ~ 2. If is however, imporlanl fo be able fo communicate 
fhe deleclion promplly, so lhal an immediate campaign fo localize fhe evenl and Irigger higher 
frequency insfrumenfs may be launched wilhin a week from fhe TDE discovery. A lol of rele- 
vanl informalion will be galhered from TDE discoveries fhanks fo fhe synergy among fulure radio, 
oplical and X-ray surveys, allowing fo: 

• make fhe eslimale of fhe jel produclion efficiency more reliable 

• differenliale differenl Iheories of jels produclion onto supermassive black holes, e.g., spin- 
powered versus accrelion-disk powered 

• discover quiescenl SMBHs in fhe dislanl Universe 

• probe fhe BH mass funclion al low masses. 
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